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ABSTRACT 

Soil moisture content can be estimated by evaluating the velocity at which sound waves 
travel through a known volume of solid material. This research involved the development of three 
soil algorithms relating the moisture content to the velocity at which sound waves moved through 
dry and moist media. Pressure and shear wave propagation equations were used in conjunction 
with soil property descriptions to derive algorithms appropriate for describing the effects of 
moisture content variation on the velocity of sound waves in soils with and without complete soil 
pore water volumes, An elementary algorithm was used to estimate soil moisture contents ranging 
from 0.08 g/g to 0.5 g/g from sound wave velocities ranging from 526 m/s to 664 m/s. Secondary 
algorithms were also used to estimate soil moisture content from sound wave velocities through 
soils with pores that were filled predominantly with air or water. 

INTRODUCTION 

Ultrasound methods involve the mechanical vibration and propagation of waves above 
about 20,000 cycles/s through various materials (Dull et al., 1964), The velocity of these 
propagating waves is affected by the nature of the material through which it is passing. Curtis 
(1982), Szilard (1982), and Kinsler et al. (1982) describe several equations and input parameters 
used in estimating the velocity of sound waves traveling through different media. Sound waves 
travel through air and water at velocities of approximately 330 m/s and 1660 m/s, respectively. 
These waves travel through solid materials such as aluminum and steel at velocities of 6,300 m/s 
and 6,100 m/s, respectively. They also travel through porous materials such as concrete, ice, and 
cork at velocities of 3,100 m/s, 3,200 rids, and 500 m/s, respectively (Kinsler et al., 1982). 
Kinsler et al. (1982) also reported a number of other solid material properties associated with 

sound wave velocties including density (kg/m 3 ), Young’s Modulus (Pa), Shear Modulus (Pa), 
Adiabatic Bulk Modulus (Pa), Poisson’s Ratio, and Characteristic Impedance (Pa* s/m). 

In-situ soils have some properties that are similar and different than those of many materials 
classically evaluated using ultrasonic sound wave velocity measurements. Soils are quite 
heterogeneous and the classically evaluated solid materials are homogeneous (Brady, 1990; 
Szilard, 1982). Soils can be compacted to a greater degree than materials such as steel and iron 
and have porosities around the range of materials such as cork and oak wood. Moisture content 
variation in soils occurs differently than in solid materials such as steel and iron, or glass and 
quartz. Texture, organic matter content, bulk density and porosity are major contributors to the 
soil water holding capacity (Brady, 1990). Additional contributions may be made by the presence 
of soil microbes, rocks, and other chemical and physical parameters. The purpose of this study is 
to perform initial investigations of the derivation of elementary and secondary ultrasound 
algorithms useful in evaluating soil moisture content variation. The study will attempt to 
incorporate soil moisture content, porosity, organic matter content, bulk density, and texture into 
existing equations that are used to estimate the velocity at which ultrasonic waves propagate 
through soils with pores that are filled predominantly with air or water. 

ALGORITHMS 

Sound waves moving through media with no boundary effects can be described by the 
sound velocity, c, which can be described by shear and pressure waves (Szilard, 1982). 
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and 
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where velocity is determined by the density, p, the modulus or elasticity, E, and Poisson’s ratio, 

G, and G is the modulus of rigidity (also called the shear modulus or torsional modulus). The 
density of a given volume of solid material is described similarly to the density of a given volume 

of soil, Pr; density of soils is equal to the mass of the dry soil, M s , divided by the entire or bulk 
volume occupied by the material as shown in Equation (3) 


ps = 


Ms 

Vtotal 


(3) 


where M s is the soil mass and Vtotal is the total volume of a given soil sample. The total 
velocity, cjoTALj of sound waves through a medium may be described as 


CTOTAL = Cpress + Cshear (4) 

Therefore, 
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Equation (5) suggests that the velocity of sound waves moving through solid media decreases with 
positive increases in p. 

Soil moisture content, 0, is defined as 
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Substituting Equation ( 6 ) into Equation (3) results in 
Mw 
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Substituting Equation (7) into Equation (5) results in 
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CTOTAL can be further defined as 
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where 


a = Vtotal 
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Equation (9) can be substituted into Equation (8) and results in 
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Equations (8) and (10) are the elementary algorithms (Eland Eli) used for estimating mass and 
volumetric soil moisture content, respectively, from sound wave velocity through soils. The 
effects of porosity on cjOTAL can be described using a modification of either of the elementary 
algorithms as shown in Equation (11) and ( 12). 
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where VsOLID is the volume of solid portion of the soil system. Organic matter and texture effects 
can also be evaluated through the measurement and evaluation of the volume of soil organic and 
mineral fractions within VsOLID- 

Szilard ( 1982) further describes the velocity of ultrasound waves through air and water as 


and 
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where cfluid is the speed of sound waves through air, y is the specific heat ratio, P 0 is the static 
pressure, Po is the static density, CpLUID is the speed of sound waves through water, k is the bulk 

stiffness modulus, and p is the density. Soil particles are considered to be in air until the water in 
the pore space flows or has similar properties as a complete pore volume of water. Equation (7) 
can therefore be substituted into Equations (13) and (14) to result in 
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Cfluid = ■ 
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Equations (15) and (16) are the secondary algorithms (SI and SII) used for evaluating c through 
porous media in air and water, respectively. Cair and Cfluid can be further defined as 


CaIR = V OmA 


(15a) 
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RESULTS AND DISCUSSION 


The El elementary algorithm (Equations 8a, b, and c) for estimating sound wave velocity 
from mass soil moisture content was used for the algorithm simulation. lOOgmof a well 
granulated soil and the preliminary data from Choi et al. (1996) were used as input data. The value 
of a was assigned a value of 0.3 and the shear propagation through the media was assumed to be 
negligible. The mass moisture content for the simulation ranged from 0.08 (wilting point) to 0.5 
(saturation) and the numerical value for the modulus of elasticity, E, ranged from 9.25E5 (moist) 

to 1.47E6 (dry). Figure 1 shows the results of plotting 0 m versus cxoTAL f° r the algorithm 
simulation. The velocity of the sound waves decreased from 664 m/s to 527 m/s with a slope of - 
312 as the soil moisture content increased from 0.08 g/g to 0.5 g/g. 

Soil moisture contents ranging from 0.08 g/g (wilting point) to 0.5 g/g (saturation) were 
estimated from the velocity of sound waves moving through the soil media using the secondary 
algorithms SI and SII (Equations 15a and b and 16a and b). The results of plotting mass moisture 

content, 0 m , versus sound wave velocity in soils with little water in the soil pores, Cair, an d in 
soils with predominantly water-filled pores, CpLUID, are shown in Figure 2. The SI algorithm 
estimates soil moisture content ranging approximately zero to 0.3 g/g and the SII algorithm 
estimates soil moisture content ranging from approximate y 0.3 to that of saturation. 

CONCLUSION 


The ultrasound algorithms developed for estimating soil moisture content from sound wave 
velocities are viable tools which can be used in this regard. The Eland EH elementary algorithms 
are capable of estimating soil mass and volumetric moisture contents as well as other soil properties 
such as porosity, organic matter, and texture. The SI and SII secondary algorithms facilitate the 
estimation of soil moisture content from sound wave velocities at moisture contents ranging from 
approximate y zero to that of saturation. The major considerations and limitations associated with 
the algorithms include the fact that the equations described by Szilard were developed for solid, 
homogeneous materials as opposed to soils which are porous, heterogeneous materials. Material 
lattice description and consistency and the modulus of elasticity are parameter considerations that 
vary greatly from one material to the next. The described algorithms suggest a decrease in the 
porous material elasticity with increased moisture content. 
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Figure 1. Mass moisture content versus soundwave 
velocity using the elementary algorithm, El. 
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Figure 2. Soundwave velocity through soils containing 
predominantly air and water using the secondary 
algorithms, SI and Sll. 
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